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Non-Reactive  Shear-Coaxial  Jets  with  and 
without  Transverse  Acoustic  Forcing 


Sophonias  Teshome 


University  of  California,  Los  Angeles 
Department  of  Mechanical  and  Aerospace  Engineering 


Motivation 


Feedback  cycle  between  liquid  rocket  engine  (LRE)  combustion  chamber  pressure 
perturbations  and  unsteady  combustion1’2 

Large  amplitude  fluctuations  in  pressure  and  combustion  heat  release  rates  => 
combustion  instability 


Flow  Heat  Release  Pressure 

Perturbation  Oscillation  Perturbation 


t 


1  Harrje,  D.T.,  and  Reardon,  F.H..  Scientific  and  Technical  Information  Office,  National  Aeronautics  and  Space 
Administration,  NASA  SP- 194,  1972. 

2Schadow,  K.C.,  Gutmark,  E.,  Parr,  T.P,  Parr,  D.M.,  Wilson,  K.J.,  and  Crump,  J.H..  19th  AIAA  Fluid 
Dynamics,  Plasma  Dynamics  and  Lasers  Conference,  AIAA  1987-1326 
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Objective 
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•  Impose  external  acoustic  perturbations,  and  examine  the  response  and  stability 
characteristic  of  shear-coaxial  injector  flow  to  pressure  perturbation 


Acoustic/Pressure 

Perturbation 


Shear-Coaxial 
Injector  Flow 


Investigate  influence  of  injector  geometry  on  flow  response  to  external  pressure 
perturbation 


Vary  the  outer-to-inner  jet  momentum  flux  ratio,  J,  under  subcritical  chamber  pressure 
condition,  i.e.,  reduced  pressure  Pr  =  0.44,  1.05 


J  = 


PoK 

Piii] 


Pr  = 


P 


chamber 


P 

critical ,  N  2 


Perineal, n2  =3. 4 MPa  (493 psi) 


•  Characterize  mixing  using  dark-core  length  measurements 


Apply  proper  orthogonal  decomposition  of  high-speed  image  pixel  intensity 
fluctuations  to  extract  spatial  and  temporal  characteristics  of  prevalent  coherent  flow 
structures 
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Schematic  of  Experimental  Facility 
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Image  of  Experimental  Facility 


Coaxial  Injector 


Piezo-Siren 


Thermocouple  Probe 


Differential  Pressure 
Transducers 
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Waveguide 


High-Speed 

Camera 


Inner 

Chamber 


Injector  Assembly 


IJ  Injector  Tube 


* - OJ  Injector  Tube 


IJ  Inlet  Port 

i 

IJ  Plenum 


PT/TC 

Ports 


OJ  Plenum 


I  Inlet  Port - ► 


* - OJ  Inlet  Port 


PT/TC 

Port 


PT/TC 

Port 
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Injector  Configuration 
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Two  types  of  outer-to-inner  jet  cross-sectional  area  ratios 
—  Large  Area  Ratio,  thin  post  (LAR-thin) 


O 


Injector 

D, 

D, 

D4 

t/D1 

AJA, 

LAR-thin 

0.70 

0.89 

2.44 

3.94 

0.13 

10.6 

SAR-thick 

1.47 

3.96 

4.70 

6.35 

0.84 

2.90 
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Waveguides 
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Acoustic  Source  - 
fits  in  circular  end 
of  waveguide 


Rectangular  enc  - 
which  connects  to 
inner  chamber 


Standing  wave 
pressure  anti-node  case) 

Waveguides  viewed  the  front  -  cross  section 


Waveguide  viewed  from  side 
(going  from  round  to  rectangular 
cross  section) 


Inner  Chamber 


Circular  end 


Acoustic  Field  Set-Up:  Pressure  Antinode 
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Pressure  antinode  (PAN)  -  condition  of  maximum  pressure  perturbation  in  the 
acoustic  field 

Piezo-sirens  forced  in-phase 

Superposition  of  quasi- ID  acoustic  waves  traveling  in  opposite  directions  =^>  PAN  at 
the  jet  location  (geometric  center  of  test  section) 


BASELINE  FLOW 

IMAGES  AND  DARK-CORE  LENGTHS 
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Baseline  Flows:  LAR-thin  Injector 


P,  =  0.44 


Baseline  Flows:  SAR-thick  Injector 


P,  =  0.44 


Pr  =  1.05 
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Dark-Core  Length  (DCL)  Measurement 


•  First  raw  grayscale  images  were  converted  to  binary  images 

•  A  contour  was  drawn  around  the  “dark-column”  in  the  binary  image 

•  Axial  length  of  the  dark-column  measured  and  defined  as  the  Dark-Core  Length,  L 


Grayscale  (8-bit)  Binary  Image  Raw  Image  w/ 
Raw  Image  Contour 


Baseline  Normalized  DCL,  LB!D1 


25 


20 


15 


-i 

10 

5 

0 


O  LAR-thiu.  Pr  =  0.44 
*  LAR-thin,  Pr  =  1 ,05 


J _ I _ I _ L 


0  5  10  15  20  25 


/ 
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Geometric  Configurations 
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LAR-thin 


SAR-thin 


LAR-thick 

SAR-thin 

SAR-thick 

LAR-thin 


Similar  I.D. 


0 


Similar  lip 
thickne  s  s/I .  D 


Similar  I.D. 


All  dimensions  in  mm 


Di 

d2 

d3 

d4 

tlDl 

AM 

0.51 

1.59 

2.42 

3.18 

1.05 

12.9 

1.40 

1.65 

2.44  3.94  0.09 

1.6 

1.47 

3.96 

4.70  6.35  0.84 

2.9 

0.70 

0.89  2.44  3.94  0.13 

10.6 

© 


L  SAR-thick 


Two  recesses 


SAR,  LAR  — >  Small,  Large  Area  Ratio 
Thick,  Thin  — »  Post  lip  thickness 
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Baseline  Normalized  DCL,  LB!D1  -  Pr  =  0.44 
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Baseline  Normalized  DCL,  LB!D1  -  Pr  =  1.05 
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Baseline  Normalized  DCL,  LB!D1  Collapse 


Pr  =  0.44 


Pr  =  1.05 


ACOUSTICALLY  FORCED  FLOW 
IMAGES  AND  DARK-CORE  LENGTHS 
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Acoustically  Forced  Flows:  LAR-thin  -  Pr  —  0.44 
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Baseline  PAN  PN 


J  =  0.1,  fF  =  3.13  kHz 


Baseline  PAN  PN 

P'Pk-Pk/Pc  =  0M%  1-39%  0.13%  P'pk- 


Baseline  PAN  PN 

JPC  =  0.01%  0.96%  0.25% 


J  =  0.5,  fF  =  3.14  kHz 


Baseline  PAN  PN 

JPC  =  0.01%  1-42%  0.23% 
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Acoustically  Forced  Flows:  LAR-thin  -  Pr  —  0.44 


J  =  11,  fF  =  3.10  kHz 


Baseline  PAN  PN 


P'pk-JPc  =  0.01%  1.04% 


0.30% 


Baseline 

=  o.oi% 


J  =  20,  Sf  =  3.11  kHz 
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Acoustically  Forced  Flows:  LAR-thin  -  Pr  —  1.05 
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Baseline  PAN 


Baseline  PAN 


PN 


P'pk-PJPc  =  0.01%  0.77% 


Ppk-JPc  =  0.00%  0.98% 


0.11% 


J  =  0.1,  Jf  —  3.16  kHz  J  =  0.5,  fF  =  3.10  kHz 


Baseline  PAN 

p'pk-pkiPc  0.01%  0.75% 


J  =  8.5,  fF  =  3.10  kHz 


Baseline 

Ppk-pJP.  =  0-01% 


PAN 

1.22% 


J  =  1.9,  fF  =  3.10  kHz 


Baseline  PAN 


PN 


Baseline  PAN 


P’Pk-PJPc  =  0.01%  0.54% 


0.05%  Ppk-pk/Pc  =  0.01%  1.02% 
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PAN  Forcing:  40  jus  Interval  Image  Sequence 
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LAR-thin,  Pr  =  0.44,  /=  0.5 
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PAN  Forcing:  40  jus  Interval  Image  Sequence 
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LAR-thin,  Pr  =  0.44,  /=  11 
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Acoustically  Forced  Flows:  SAR-thick  -Pr  =  0.44 


J  =  o.l,  fF  =  3.10  kHz 


Baseline  PAN  PN 

p'pk-JPc  =  0.02%  1.27%  0.19% 


J  =  0.5,  fp  =  3.04  kHz 


Baseline  PAN  PN 

Ppk-pk/Pc  =  0.02%  1.42%  0.28% 


J  =  2.1,  fF  =  3.07  kHz 


Baseline  PAN  PN 

P'pk-pk/Pc  =  0.02%  1.06%  0.14% 


J  =  5.7,  fF  =  3.11  kHz 
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Acoustically  Forced  Flows:  SAR-thick  -Pr  =  0.44 


Baseline  PAN 

■p'Pk-JPc=  0.02%  1.11% 


PN 

0.23% 


Baseline  PAN 

P'pk-pJPc  =  0.02%  1.48% 


J  =  15,  fp  =  3.04  kHz 


Baseline 

?U-,  k/P*  =  002% 


J  =  21,  fp  =  3.11  kHz 
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Acoustically  Forced  Flows:  SAR-thick  -Pr=  1.05 


J  =  0.1,  /f  =  3.12  kHz 


Baseline  PAN  PN 


Ppk-JPc  =  0-01%  0.67%  0.40% 


J  =  0.5,  fF  =  3.00  kHz 


Baseline 

P'pk-JPc  =  0.01% 


J  =  2.1,  fF  =  3.04  kHz 


J  =  5.2,  /f  =  3.08  kHz 
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Acoustically  Forced  Flows:  SAR-thick  -Pr=  1.05 


Baseline 


PAN  PN 


p'pk-JPc  =  0.01%  0.99% 


0.19% 


Baseline  PAN 

PU-vk/Pc  =  0.01%  1.22% 


PN 

0.25% 


J  =  14.  fF  =  3.04  kHz 


Baseline 


PAN  PN 


Ppk-pk/Pc  =  0.°!%  1.11% 


0.13% 
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PAN  Forcing:  Normalized  DCL,  LPAN!LB 
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PN  Forcing:  Normalized  DCL,  LPN!LB 
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APPLICATION  OF  POD  ON 
IMAGE  PIXEL  INTENSITIES 
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Proper  Orthogonal  Decomposition 
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•  Proper  Orthogonal  Decomposition  (POD)  or  Principal  Component  Analysis 
(PCA)  was  used  for  extracting  dominant  dynamical  processes  embedded  in  high¬ 
speed  images. 

•  A  time-resolved  set  of  images  A(x,t)  can  be  represented  as  a  linear  combination 
of  orthonormal  basis  functions  (j)k  (aka  proper  orthogonal  modes)12 : 

M 

A(x,t)  =  ^cik(t)</>k(x) 

k= 1 

where  ak(t)  are  time  dependent  orthonormal  amplitude  coefficients  and  M  is  the 
number  of  modes 


•  Main  idea:  POD  modal  amplitudes  capture  the  maximum  possible  “energy”  in 
an  average  sense3,  i.e., 

X(a,(0«,(0)>Z(W*(0) 

k  k 

where  bk(t)  are  the  temporal  coefficients  of  a  decomposition  with  respect  to  an 
arbitrary  orthonormal  basis  y/k. 


1  Chatterjee,  A.  Current  Science,  Vol.  78,  No.  7  (2000) 

2  Arienti,  M,  and  Soteriou,  M.C..  Phys.  Fluids  21,  1 12104  (2009) 

3  Berkooz,  G.,  Holmes,  P.,  and  Lumley,  J.L..  Annu.  Rev.  Fluid  Meek  25.  539  (1993) 
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Construction  of  Data  Set 


iA*Q  I  eii^^vin  .  fiiivli O*8  i 


•  First,  form  a  row  vector  consisting  of  all  pixel  intensity  values  of  each 
snapshot  image  (with  resolution  of  n  rows  by  m  columns)  in  order  of  increasing 
columns,  then  increasing  rows 


Image  Frame 


columns 

•  Then,  combine  all  such  row 
vectors  for  N  sequences  of  image 
frames  resulting  in  a  matrix  A 
consisting  of  N  rows  by  P  =  n  x  m 
columns  of  intensity  values. 


\7 

t  ^ 

P  =  n  x  m 

-  - ► 

pixel  intensities 

TV  time 
steps 


J 
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Orthogonal  Decomposition  Technique 


•  Eigenvalue  decomposition  or  singular  value  decomposition  (SVD)  can  be  used 

•  SVD  preferred  since 

1 .  Applicable  to  non-square  matrices  (most  likely  the  case) 

2.  Decomposition  matrices  are  orthogonal 

3.  Subroutine  readily  available  in  MATLAB® 

•  Subtracted  temporal  mean  of  A  =>  matrix  of  intensity  fluctuations  A 


•  Applied  SVD 


W  r-p  rp 

=>  A  =  USVT  =  QVT 

Columns  of  Q  ~  ak{t)  contain  temporal  information 
Columns  of  V  ~  (f)k(x)  contain  spatial  information 


N*N  « 

Orthogonal  Matrix  of  Left 
Singular  Column  Vectors  of  A 


NxP 

Orthogonal  Matrix  of  Right  Singular 
Column  Vectors  of  A  <=>  proper 
orthogonal  modes  (POM) 


Nx  N 

Diagonal  Matrix  of 
Singular  Values 
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Results  -  Baseline  Low  J  Flow  at  Pr  =  0.44 


LAR-thin  ,  Pr  =  0.44,  J=  0.5 

Snapshot  Average 


Amplitude  information 
contained  in  singular  values 


Antisymmetric  Structures 

Identified  with  Characteristic 
Frequencies 


10  10 
Mode 


Power  Spectral 
Densities  (PSD) 
of  Temporal 
Coefficients  of 
POMs  1  and  2 
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Results  -  PAN  Low  J  Flow  at  Pr  =  0.44 


•  LAR-thin,  Pr  =  0.44,  J  =  0.5,  forcing  Frequency,^  =  3.14  kHz 


Snapshot 


Symmetric  Structures 

Identified  with  Characteristic 
Frequency  at  fF 


Results  -  PN  Low  J  Flow  at  Pr  =  0.44 


•  LAR-thin  ,  Pr  =  0.44,  J  =  0.5,  forcing  Frequency,  fF  =  3.14  kHz 


Antisymmetric  Structures 

Identified  with  Characteristic 
Frequency  at  fF 


/(kHz)  /(kHz) 
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Cross-Power  Spectral  Density  (CPSD) 


CPSD  yields  the  FFT  of  the  cross-correlation  of  the  temporal  coefficients 

cov(a.a) 

CPSD  =  y  "~,mk 


k~° 

Magnitude  and  phase  plots  used  to  determine  existence  of  propagating  structures 


LAR,  Pr 


0.44,  J=  0.5 


Baseline 


PAN 


Confirmed  Downstream 
Propagating  Flow 
Structures 


/(kHz) 
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Sample  Animation  -  PAN  (fF  =  3.14  kHz) 


•  LAR/V  =  0.44,7=0.5 


Raw  images  ROMs  1  and  2 
Time  -  0.00  ms  Time  =  0.00  ms 


POM  1  POM  2 

Time  =  0.00  ms  Time  =  0.00  ms 


t 

Superposition  of  POMs  1  and  2  Resulted  in  Downstream  Propagating  Structures 
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CPSD  Mi-ifHiimlt 


Baseline  POM  &  CPSD:  LAR-thin  -Pr  =  0.44 


Baseline  POM  &  CPSD:  LAR-thin  -Pr  =  1.05 


/=  0.5 


fik Hi> 


1.9 


/(kHz) 


5.0 


/(hBz) 
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N  on-Dimensional 
Spatial  Growth  Rate 


Previous  Works  on  Jet  Instability 


Michalke  and  Hermann  (1982)  did  linear,  inviscid  instability  analysis  of  a  circular  jet 


with  coflow 

Showed  that  with  increasing  coflow  velocity, 

—  Helical  disturbances  more  unstable  than  axisymmetric  ones  farther  downstream  of  exit 

—  Jet  flow  becomes  less  unstable,  but  spectrum  of  spatial  growth  rate  becomes  broader  and  the  peak  shifts 

to  higher  frequencies  UU1„ 


Helical 


Axisymmetric 


>  Non-Dimensional  Frequency 


•  Dahm  et  al.  (1992),  Wicker  and  Eaton  (1994)  conducted  experimental  investigation  of 
large-scale  vortex  structures  in  the  near  field  of  coaxial  jets 

—  For  outer-to-inner  jet  velocity  ratios  greater  than  one,  found  that  coherent  structures  in  the  outer  shear 
layer  dominate  those  in  the  inner  shear  layer 
—  At  large  axial  distances,  shear-layer  vortices  exhibit  helical  structures 
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Dahm  et  al ,  JFM  1992 


{’PSI>  Miiiinilmlf 


Baseline  POM  &  CPSD:  SAR-thick  -Pr  =  0.44 


CPSli  IViHgnitiKle 


Baseline  POM  &  CPSD:  SAR-thick  -Pr=  1.05 


bfjici  «  .  Sift¥lieP  i  Hurt*  t 


Inner  Shear-Layer  Structure  Velocity  -  Pr  =  0.44 


LAR-thin  SAR-thick 
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Inner  Shear-Layer  Structure  Velocity  -  Pr  =  1.05 
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Conclusions 
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o  Baseline  flow  normalized  dark-core  lengths,  LBIDX\ 

•  Power-law  curve-fit  of  LB/Dl  variation  with  J  showed  that  Pr  =  1 .05  cases  had  a 
stronger  dependence  (J° 50  for  the  LAR-thin  and  54  for  the  SAR-thick);  Pr  =  0.44 
cases  varied  as  J~° 39  for  the  LAR-thin  and  35  for  the  SAR-thick. 

•  For  a  given  J  and  same  Pr,  SAR-thick  flows  had  smaller  LB/DX  than  LAR-thin  flows 

•  Considering  the  LAR-thick  and  SAR-thin  injector  flow  data  at  both  Pr,  it  can  be 
deduced  that  LB/Dl  increased  with  AJAi  while  it  decreased  with  J  and  t/D{ 

=>  injector  geometry  also  plays  an  important  role  in  addition  to  J 

o  Acoustically  forced  flow  normalized  dark-core  lengths,  LPANiILb  or  LPN/ILB. 

•  LAR-thin  injector 

■  Both  LPANiILb  and  LPN/ILB  less  than  one  and  approached  one  with  increasing  J 
=>  The  flow  became  less  sensitive  to  forcing  with  increasing  J. 

•  SAR-thick  injector 

■  For  low  J  (<  10)  flows,  LPAN/ILB  stayed  around  one 

■  Lpn/ILb  approached  one  with  increasing  J 

=>  Flow  in  recirculation  zone  only  effective  in  dampening  the  effect  of  PAN 
forcing  at  lower  J 
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Conclusions  (cont’d) 


o  Dynamic  characteristics  of  baseline  flows: 

•  LAR-thin  injector: 

■  inner  shear  layer  formed  immediately  downstream  of  the  exit 

■  baseline  flow  inner  shear  layer  characterized  by  antisymmetric  flow 
structures  that  indicated  the  presence  of  helical  instabilities 

■  peak  frequencies  of  the  magnitude  spectra  became  broader  and  shifted 
to  higher  frequencies  with  increasing  J 

•  SAR-thick  injector: 

■  region  immediately  downstream  of  the  thick  inner  jet  post  formed  a 
flow  recirculation  zone  delaying  formation  of  shear  layer  between  the 
bulk  inner  and  outer  jets 

■  peak  frequencies  stayed  at  low  frequencies  despite  increasing  J 
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Conclusions  (cont’d) 


■  S>iv1i Q*  I  eimvli 


o  Dynamic  characteristics  of  forced  flows: 

•  LAR-thin  injector: 

■  PAN  forcing  produced  symmetric  flow  structures  in  the  inner  shear  layer  region 
for  low  J  (<  5)  flows  at  both  Pr 

■  Magnitude  spectra  showed  dominant  peaks  identical  to  the  forcing  frequency 

■  Higher  J  (>  5)  flows  under  PAN  forcing  produced  antisymmetric  structures  that 
resembled  that  of  baseline 

■  Magnitude  spectra  retained  baseline  flow  spectral  characteristics 

■  peak  frequencies  of  the  magnitude  spectra  became  broader  and  shifted  to  higher 
frequencies  with  increasing  J 

■  PN  forcing  produced  sinusoid  disturbance 


SAR-thick  injector: 

■  PAN  forcing  produced  no  appreciable  response  in  the  inner  shear  layer  of  the 
low  J  (<  5)  flows,  at  both  Pr 

■  Spectra  showed  low  frequency  peaks,  and  no  peaks  at  forcing  frequencies  were 
detected 

■  at  higher  J  (>5)  and  both  Pr,  PAN  forcing  produced  symmetric  structures  below 
the  recirculation  zone  with  peaks  in  the  spectra  at  the  forcing  frequencies 

■  PN  forcing  also  produced  appreciable  response  at  higher  J 
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Acoustic  Field  Set-Up:  Pressure  Node 
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Pressure  node  (PN)  -  condition  of  minimum  pressure  perturbation  in  the  acoustic  field 
Piezo-sirens  forced  out-of-phase 

Superposition  of  quasi- ID  acoustic  waves  traveling  in  opposite  directions  =^>  PN  at  the 
jet  location  (geometric  center  of  test  section) 
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